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In the framework of Relativistic MOND theory (TeVeS), we show that a late time back- 
ground ACDM cosmology can be attained by choosing a specific F(fi) that also meets 
the requirement for the existence of Newtonian and MOND limits. We investigate the 
dynamics of the scalar field <f> under our chosen F(fi) and show that the "slow roll" 
regime of <j> corresponds to a dynamical attractor, where the whole system reduces to 
ACDM cosmology. 
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1. Introduction 

Current as tronom ical observations on cosmological scales (CMB anisotropy 
supernovae ISEEl and SDS^) reveal that our universe is spatially flat, with about 
two thirds of the energy content resulting from what is referred to as dark energy. 
This energy has a negative pressure to account for the accelerated expansion of 
the universe. While, on smaller, galactic scales, the rotational curves of galaxies 
strongly indicate that the biggest contributions to their mass density arises from 
non-luminous matter, which has given rise to speculations on the existence of dark 
matter^. 

On the other hand, it is also interesting to inquire if we can solve the current 
observational puzzles, both dark matter and dark energy, by modifying Einstein 
gravity in spite of its many successes, in particular, the solar system tests. MOND 
theory is a striking modification of the Newtonian laws of motion and explains the 
rotational curves with amazing accuracy without introducing any dark matter. 
After some initial problems, a consistent relativistic extension of the MOND theory 
has recently been proposed by BekensteirJ^, which has spurred extensive studies on 
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its successes and pitfalls ^ l ^ l ^ | 13 | 14 | 15 | 16 | 17 | 18 | j n that extension, general relativ- 
ity has been modified by including two scalar fields (one dynamical while the other 
is non-dynamical), one vector field. In addition, the theory contains one arbitrary 
function whose form is dictated in part by its primary purpose, i.e., to address the 
data on the rotational curves of galaxies, or equivalently, the dark matter problem. 
Consistency with data from the Bullet cluster indicates that eV scale hot neu- 
trinos must also be included as explicit dark matter candidate^! The remarkable 
features of this theory are that it can explain the galaxy rotational curves without 
introducing dark matter (just as its non-relativistic counterpart, MOND theory) 
and at the same time reduce to general relativity in the appropriate limits. The 
price one pays is the additional complication of the whole system: two scalar fields, 
one vector field and gravitational field. 

A natural question to ask is if the new framework also works for cosmology? Can 
it help to tackle the dark energy problem? We will show in this paper that under 
appro priate condition s, the scalar field could mimic the behavior of the dark energy 
models ^ | 22 | 23 | 24 | 25 | Specifically, in the "slow roll" regime, it will be seen to reduce 
to the standard ACDM model. Although the current formulation of TeVeS is still 
not as mature as the dark matter models, it has the potential, to naturally solve 
the cosmic puzzles currently ex plai ned by dark energy and dark matter without 
invoking either of them explicitly^ 1 . 



2. Basic setup 

We begin with a brief outline of Gravity in the TeVeS framework. Following 
BekensteirP, the action for the whole system could be written as 



where 



S Sg + S s + S v + S m (1) 

S g = (WirG)- 1 J g a0 R aP V=gd 4 x (2) 

S s = -\J [ct 2 /i q/3 0, q </>,/3 + ^Gr 2 c7 4 F(kGa 2 )]y/^d 4 x (3) 

Sv = -3^3 J [rt^AH - 2(\/K)(g^ii v + 1)] V^9d 4 x, (4) 



s m = j c($ MV ,r,ffr,-~W-9<r* (5) 

are the actions for gravity, scalar field, vector field and luminous matter respectively. 
In the framework of TeVeS, there are two different frames, physical frame corre- 
sponding to the "real universe" and Einstein frame. The metric in physical frame 
(ffjui/) relates that in Einstein frame (g^v) by g a p = e^ 2 ^(g a p + ilailp) — e^t&aiip. 
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Here il Q is a timelike 4- vector field, cj> and a are respectively the dynamical and the 
non-dynamical scalar helds, F is a free function to be specified by dynamics while 
k and K are two positive dimensionless parameters. Varying the action, one can 
arrive at the following equations of motion for the vector, scalar and gravitational 
fields^ 

K (^ a ^l, +il a il y ii^. \ + 8vrGa 2 [i^0 i/3 <? Q >, 7 +U°(iC^) 2 ] 

= 8ttG(1 - e~ i(i> ) [gfUPftf + H^WT^p] (6) 



- /xF( M ) - -Li 2 F'( f i) = k£ 2 h af3 ^i3 (7) 
[H (kfh» v <l>^ v ) h a ^, a ]. fi = kG[g al} + (1 + e-^)^]f aP . (8) 



G a p = 8ttG 



T afJ + (1 - e- A *)WT Ka il m + r af3 



e 



a/3 



(9) 



where 



T a fJ = & 
1 



Gr 2 a 4 F{kGa z )g a0 



(10) 



e af3 = K^Ufa^n - ^g° T g» v il M tt [TM g aP ) - Xilattp (11) 

T^u is the energy momentum tensor of ordinary matter in the physical coordinate 
system. h af> = g afi - iPif 3 and fi = kGa 2 . 



3. Cosmological dynamics in Einstein frame 

In a cosmological setting, the symmetries of the FRW universe will simplify the 
above equations considerably. It is worth recalling here an underlying assumption 
in such applications. One interprets the fields in the problem, say the scalar field, 
as consisting of both a spatially homogenous (time dependent) part and an inho- 
mogeneous one. On cosmological scales the spatially inhomogeneous part may be 
neglected in the first approximation while at galactic scales, the inhomogeneous 
part plays a prominent role and the homogenous part is negligible (quasi-static 
approximation). Henceforth, we will restrict ourselves to cosmological scales, more 
specifically to the flat universe case, and consider the line element: 



ds 2 = -dP + a(t) 2 dx 2 



(12) 
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where dt = e^dt and a = e~^a. It is clear that the physical frame and Einstein 
frame become the same once the is constant 1. The vector field in this case could 
be chosen as il M = 5!£ and the scalar field depends only on time t. In Einstein frame, 
the equations of motion and the field equations have simpler forms though their 
dynamical evolution are essentially the same as those in physical frame (they are 
related by a non-singular transformation). So, in this section, we will analyze the 
dynamics of the system in Einstein frame. 



d\ 2 8ttG 



a J 3 
4ttG 



(pM + Ps) (13) 



(pm + P4> + 3pm + 3p<t,) (14) 



a 3 

where pu = pe~ 2 ^, Pm = pe~ 2 ^ are respectively the density and pressure of non- 
relativistic matter, p^ and p^ are the effective energy density and pressure of the 
scalar field <fi, defined as: 

^ kG + U 2 k 2 G W U 2 k 2 G [ ' 

= pft_ p? p 2 [F(p) + pF'(p)} 

P0 kG U 2 k 2 G lMJ At 2 k 2 G 

where p relates to <p via the a equation (EqjTJ): 

pF(p) + l -p 2 F'(p) = 2k£ 2 <p 2 (16) 

Here we need to emphasize, as suggested in 03 th a t the above definitions are held 
in the effective sense due to the fact that pM and pm terms contain <fi. But in the 
"slow roll" regime (<j) ~ as we will see later), they are well defined as the <j> is 
merely a constant factor. The equation of motion of the scalar field 4> in a FRW 
background is given by 

4>+(3H+^)<p+—{p M + 3p M )=0 (17) 
p Zp 

where H = a/a and dot denotes the derivative with respect to t. It is worth nothing 
that, under the specific form oiF(p) we choose in Eq. (|2T|) . Eq. (fi~7f is consistent with 
the "slow roll" approximation (j) ~ under two different situations: First, as in the 
usual scalar field dark energy dynamical system, the matter density Pm decreases 
to zero at late times, which leads to a dominated universe; Second, - « 1, which 
lead to a universe with both pM and p§. In next section, we consider the second 
possibility when fitting the model to the Supernovae data. On the other hand, if we 
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assume that <fi varies slowly, then we can consider the two energy components to be 
approximately adiabatic, and energy momentum conservation leads to 

^± = -3% l+Pl ) (18) 
at a 

where the subscript i denotes M and <j>. To get a closed system of equations, we also 
need to specify the equation of state of the energy/matter content. This equation 
of state is specified by w — p/p, which is for non-relativistic matter. While for 
the scalar field, this equation of state is given by its equation of motion, i.e..Eq. (|17p 
which can be expressed as, 



_ U 2 kp4> 2 - P 2 F (p) _ F(p) + pF'(p) 
W0 M 2 kpj> 2 + p 2 F(p) 3F(p) + pF'(p) 

In the "slow roll" regime, <j) ~ <f> ~ and ~ — 1. It is worth noting that 
from Eq. (fTO|) . the equation of state could approach —1 from either greater than 
-1 (quintessence case) or less than -1 (phantom case) if we choose the form of F(p) 
appropriately. In this paper, we focus on the case with > — 1. Eqs. (|13M14[) are 
expressed in terms of the scale factor in the Einstein picture. In the following, we 
will continue to work in Einstein frame to investigate the dynamics of the scalar 
field. 

From the above discussion, it is clear now that we will get the evolution of the 
whole dynamical system if the form of F (p) is specified. However, at present there 
are no theoretical arguments in favor of any specific choice, thereby providing a lot 
of freedom in this regard. The constraints on the form of F(p) are phenomenological 
in nature motivated by the condition that the correct physics is obtained. Besides 
the toy model F(p) suggested by Bekenstein, there are some others forms suggested 
to account for the better fit to 

data™. 

In this section, we will consider the the 
choice of the form of F(p) that will lead to an acceptable cosmology in addition to 
accommodating the MOND theory and Newton's Law. From Eq. (fTS"]) . the positive 
energy condition, p^ > and the non-vanishing of p^, as (f> — ► may restrict the 
freedom in the choice of the form of F(p). On the other hand, the existence of 
MOND and Newtonian limits also give constraints on the form of F(p) aJ^ 

d(p 2 F(p)) k 2 p 2 

dp - - 64^(1-/0"^ ) (20) 

with f(p) an arbitrary function of p with the requirement of non- vanishing /(0) 
and m — 1 is the case proposed by BekensteirM Therefore, we suggest a form of 
F{p) that meets the above requirements as 



3 p(4 + 2p- 4p 2 + p 3 )+2 m[(l - pf] i a 
8 p 2 



F (W = o ~% + ~a ( 21 ) 
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where a is a dimensionless constant. When a = 0, the above reduces to the form 
irM Then, from Eq. (fl6]) , we can obtain the relation between fi and <\> as 

^(1 + p-3p 2 + ^ - Y^~T? = AM2 ^ 2 (22) 

For a consistent cosmology, we require 2 < p < oo. Note also that = when p = 2, 
for which the energy density Eq. (fTS"]) becomes p^ = 4 ^ 2 ", 2G and the equation of state 
reduces to, w$ = — 1 at late times, which corresponds to a de Sitter phase. This 
not only provides an acceptable cosmological state but also ensures the consistency 
of our previous assumption </> ^ 0. Next, we show that such a phase corresponds 
to a dynamical de Sitter attractoJ^. To do this, we introduce the dimensionless 
variables x = <fi, z = to4> and N = In a with to a scaling constant of dimension t. We 
can rewrite the Eqs. (113lll8[) in terms of x, z and N and linearize them around the 
critical point (x, z) = (x, 0), then, we arrive at the following system of equations 



^ = Z (23) 

dN t 0y /2Tra/3k 2 £ 2 

dN 

It is easy to see that the eigenvalues of the coefficients matrix of Eq.(|23p are (—3, 0) 
which indicates that the critical point is stable, i.e. a dynamical attractor. One 
comes to the same conclusion by noting that (j) — leads to a minimum of the 
energy density p$. In Figs. ( 1 ) and ( 2 ), we plot the numerical results for the 
dynamical system defined above. Note that our intention here is merely to illustrate 
the possibility of a consistent cosmology and not to fit the exact observational data. 
So, for convenience, we have set the parameter a = 0.01 and the rest to unity. Since 
the attractor corresponds to all <j> with <j> = 0, we can choose </> ~ so that e~ 2 ^ ~ 1, 
which is indicated by the previous discussion. In our numerical analysis, we chose 
4> = 0.00001 and the initial <f> from 0.001 to 0.005 with an increment of 0.001. We 
should point out that the above conclusion won't change even if we choose other 
values for (p. Then the physical frame and Einstein frame are different by a factor 
of e~ 2 ^. In Fig. (3), we plot the corresponding evolution of the and with 

respect to fi given the form of F(/x) in Eq. ((2"Tj) . 



4. Cosmology in physical frame 

The analysis in preceding section is carried out in the Einstein frame. But we need 
to work in physical frame if we want to compare the theory with observation data. 
The two frames are mathematically equivalent to each other although the physi- 
cal interpretations are different. With a simple transformation, one could convert 
the the equations from Einstein frame to physical frame. Under the "slow roll" 
approximation^ ~ ~ 0), the system can be casted as: 
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Fig. 1. The evolution of cosmic parameters for matter(indigo curve) and <f> field (pink curve). 
From the plots, one can find the curves with different initial conditions converges very quickly to 
a common track corresponding to the attractor. 
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Fig. 2. The evolution of the equation of state of <j> field for different initial 



(24) 



4ttG^ _ 



(25) 
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^L = -3^ +Pl ) (26) 
at a 

Where the prime denotes derivative with respect to t and we have used = —p^ 
and e -2 ^ ~ const under the condition <j> ~ 0. 

In the spirit of MOND, p is the density of baryons. Comparing Ea.(|24p to the 
standard ACDM cosmology, we can find that the term pe~ 4 ^ plays the equiva- 
lent role of non-relativistic matter (baryons and dark matter). That is to say, 
_ p — p(i _ e -40^ can ^ e considered as an "effective" dark matter den- 
sity. We introduce f2 b = p/p cri u = P^e" 2 ^ / 'p crit and Q M = pe^ / p crit as the 
energy density fraction of baryon, dark energy and non-relativistic matter respec- 
tively, where p cr u = 3H 2 /8irG is the critical energy density. From Ea.(|16p. one can 
observe that for any choice of F[fx), a is a constant as (f> ~ and therefore is 
constant in the "slow roll" regime. As a result, the above system could be reduced 
to ACDM universe with the relation 

<f>= 1 I M§ L ) (27) 
4 \l M 

From the discussion in the previous section, we know that <f> = with <f> at any value 
is a dynamical attractor. From Eq. (|27p . we can fix the cj> value via observation, i.e. 
the ratio of From Big Bang Nucleosynthesis^ 21 , we can give a constraint on 
the baryon density (p and hence Q;,), while from snieP, we can constrain the total 
density of non-relativistic matter(/5e -4< ^ and hence Qm)- Combining the two, we 
can give a constraint on <f> as <fi — —0.46 ± 0.03 and the results are shown in Fig. (4). 
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Fig. 4. Constraint on <f> from Gold Sample SNIa data. The fj, axis in this plot is the distance 
modulus. 



5. Discussion 

In this paper, we show a possible late time background ACDM like cosmology in 
the framework of TeVeS theory by choosing the F(/x) as a specific form that also 
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meets the requirements of Newtonian and MOND limits. We show the dynamical 
system has a late time de Sitter behavior and the scalar field at the de Sitter 
attractor should be —0.46 ±0.03 so as to be in agreement with current astronomical 
observations. 

On the other hand, the framework of the TeVeS theory involves an arbitrary 
function in addition to many parameters as well as auxiliary fields. This not only 
provides much freedom in matching the theory to observations, but also adds uncer- 
tainties. Clearly, they must be constrained from further phenomenological investi- 
gations of the testable predications of the theory. It would be interesting to see if the 
model can be made compatible with all the astronomical data with a unique choice 
for the arbitrary function. Are there any theoretical constraints on this function? 

The analysis in this paper is aimed at suggesting a possibility to reduce the 
complicated dynamical system of TeVeS into the customary ACDM under slow 
roll approximation and thus get some feedback from the cosmological observational 
constraints on Q^/CIm- The choice of the function F(/j,), though without any deep 
theoretical basis, phenomenologically explains both dark matter and dark energy 
in the same framework. It remains to be checked in future how this form of the 
function F(/i) can be embedded in the larger context of an appropriate function 
that explains all the data. 
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